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Spectrum analysis poses some of the most important problems in quantum mechanics. 
Success of nonrelativistic potential models of quark confinement [1] reattracted attention to 
the bound-state problem in relativistic physics as well [2,3]. Several attempts have been made 
to describe relativistic systems in a central field due to a Lorentz-vector and Lorentz-scalar 
interaction within the framework of the Dirac equation. However, for almost all potentials this 
equation is not exactly solvable which compels to resort to some approximation methods. 

A number of such approaches to solving the radial Dirac equation in analytical expressions 
have been developed, including, in particular, the use of the WKB-method [4,5], the hypervirial 
and Hellman-Feynman theorem [6,7], the 1 /iV-expansion [8-17], the algebraic approach [18], the 
method of Regge trajectories [19-21] and various perturbation schemes [21-26]. 

Despite such a variety of methods one of the most popular techniques is still remained 
logarithmic perturbation theory [27-32]. This technique involves reducing the Dirac equation, 
which becomes a pair of coupled first-order differential equations in central-field problems, to 
a nonlinear Riccati equation. In the case of ground states, the consequent expansion in a 
small parameter results in handy recursion relations which permit us to derive the high order 
corrections to energy eigenvalues and eigenfunctions. It should be emphasized that high orders 
of expansions are needed for applying modern summation procedures because the obtained 
perturbation series are typically divergent. However, when radially exited states are considered, 
the standard approach becomes extremely cumbersome and, practically, inapplicable. It is 
caused by factoring out zeros of the wave functions, which, in addition, are not the same zeros 
for small- and large-components of the Dirac spinor [33] . 

On the other hand, it is known, that the radial quantum number, n r , most conveniently 
and naturally is introduced in consideration by means of quantization conditions, as in the 
WKB-approach [34,35]. However, since the WKB-approximation is more suitable for obtaining 
energy eigenvalues in the limiting case of large quantum numbers and the perturbation theory, 
on the contrary, deals with low-lying levels, the WKB quantization conditions need change. 

Recently, a new technique based on a specific quantization condition has been proposed to 
get the perturbation series via semiclassical expansions within the one-dimensional Schrodinger 
equation [36]. For the Dirac equation, on performing the scale transformation, r — >■ h 2 r, the 
coupling constants appear in common with powers of Planck's constant, h, thus implying the 
possibility to obtain perturbation expansions in a semiclassical manner in this case, too. 

The objective of this letter is to develop the explicit semiclassical treatment of logarithmic 
perturbation theory for the bound-state problem within the framework of the Dirac equation 
and to describe a new procedure for deriving perturbation corrections through handy recursion 
formulae having the same simple form both for ground and exited states. 

The proposed technique can be regarded as a further investigation of a part assigned to a 
rule of achieving a classical limit in construction of semiclassical methods. In addition to the 
rule, h — > 0, n r — > oo, I — > oo, hn r = const, hi = const, required within the WKB-approach; 
and conditions, H — > 0, n r = const, I — > oo, Tin r — > 0, hi — const, which are applied within 
the method of 1/N-expansion [17]; here we address ourselves to the alternative possibility: 
h — > 0, n r = const, / = const, hn r —> 0, hi — > 0, that results in the explicit semiclassical 
treatment of the logarithmic perturbation theory. 

(1) Method. We study the bound state problem for a single fermion moving in an attractive 
central potential. This potential contains both the time component of a Lorentz four-vector, 
V(r), and a Lorentz-scalar term, W(r), which, in general, have a Coulomb-like behaviour at 
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the origin 

1 OO 1 oo 

V(r) = Vi r l , W(r) = , (1) 



r i=0 r i=0 

though the case V = or W = is permissible, too. In what follows, a scalar potential 
will be included in the mass term m(r) by analogy with "dynamical mass" models of quark 
confinement [2,3] 

W(r) 

m(r)=m+—-^. (2) 

cr 

Then the Dirac radial wave equations have the form 

hF'Or) - —Fir) + -\E - V(r) - m(r)c 2 ]G(r) = , 
r c 

TiG'ir-) + ^G(r) - -[£ - V(r) + m(r)c 2 ]F(r) = . (3) 
r c 

Here F(r) and G(r) are the small and large components of the wavefunction of a particle and 
E is its total energy, x — s ' + \) for j = / — |, with s = ±1 denoted the sign of \- 

Eliminating F(r) from the system (3) and performing the substitution, R{r) = hG'{r)/G{r), 
for the logarithmic derivative of large component we then arrive at the Riccati equation 

kR'(r) - hQ(r)R(r) + R 2 (r) 

= ^ h Q(r) h m z (r)c z S - V(r)J z , 

r z r &■ 

with Q(r) = [m'(r) - + m(r) - V(r)}. 

As was above pointed out, logarithmic perturbation theory does involve coupling constants 
in common with powers of Planck's constant. Therefore we attempt now to solve eq.(4) in a 
semiclassical manner. Taking into account the leading orders in h of the quantities E ~ 
%X ~ ^5 from Riccati equation (4) we have 

1 oo 

n i=0 

oo 

Q(r) = ^ 2 ^Q l (r)^, (5) 

i=0 

1 oo 
' i i=0 

Keeping in mind that relativistic mechanics must go over to nonrelativistic one as the speed 
of light tends to infinity and quantum mechanics must go over to classical one as fi — > 0, the 
correlation between these constants has to be set. By analogy with quantum electrodynamics 
the foregoing consideration will be carried out under the condition he ~ 0(1) [37]. Moreover, 
for simplicity we put he — 1. 

On substituting the expansions (5) into the Riccati equation (4) and comparing coefficients 
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of the different powers of h, one obtains the following hierarchy of equations 
Rl(r) = m 2 



Roi^R^r) = E [V(r) - E ± \ +mW(r) , 
R[(r) + 2R (r)R 2 (r) + R\{r) - R (r)Q (r) 



+ 2E 1 V(r) - El - 2E E 2 + W 2 {r) - V 2 {r) , 



(6) 



R' k ^{r) + Y, R Ar)Rk-Ar) 

3=0 

k 

= 2E k _ 1 V{r)-J2E j E k - j - 

j=0 



k—2 
3=0 

-Qk-s(r) , 



where 
Qo(r) 

Qk(r) = 



1 



En + m 



[W(r) - V\r) 



E + m 



k-i 



£ Q 3 {r)E k _ 3 + Q k ^(r) [W(r) - V(r) 

b=o 

For nodeless states this system can be solved straightforwardly. However, when radial 
excitations are described with standard technique the nodes of the wave functions need to be 
factored out first and consideration becomes extremely cumbersome. We intend to circumvent 
this difficulty by making use of the quantization condition. Its fundamental idea that stems 
from the WKB-approach [34,35] is well known as the principle of argument in the analysis of 
complex variables. Being applied to the logarithmic derivative, R(r), it means that 



j R{r) dr = hN, 



(7) 



where iV is a number of zeros inside a closed contour. 

This condition is exact and is widely used for deriving the high-order corrections to the 
WKB-approximation and the 1/iV-expansions. There is, however, one important point to note. 
The radial and orbital quantum numbers, n r and /, correspondingly, are specific quantum 
notions and need be defined before going over from quantum mechanics to classical physics. 
Therefore the quantization condition (7) must be supplemented with the rule of achieving a 
classical limit that stipulates the type of semiclassical approximation. 

In particular, within the WKB-approach the passage to the classical limit is implemented 
using the rule 

h — > 0, n r — > oo, I — > oo, hn r = const, hi = const, (8) 

whereas the 1/iV-expansion, being complementary to the WKB-method, requires the conditions 
[17] 

H — > 0, n r = const, / — > oo, hn r — > 0, M = const. (9) 

The semiclassical treatment of logarithmic perturbation theory proved to ivolve the alter- 
native possibility: 

H — > 0, n r = const, / = const, hn r — > 0, HI — > 0. (10) 

Notice that the part of this rule concerned the orbital quantum number, /, differs from one 
used within the WKB-approach and the 1/ iV-expansion method. In our consideration this part, 
implying the first order in % for the quantity hx, has been used in deriving the system (6). 
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The remaining part of the rule respects the radial quantum number. Due to (10) the right- 
hand side of the equality (7) has the first order in H and, hence, on substituting the expansion 
(5) the quantization conditions (7) takes the form 

±J Rt (r)dr = N5 hl , (11) 

where the Kronecker delta Sij is used. 

Before proceeding further, we must specify the quantity N. In contrast to the WKB- 
approach and the 1/N- method, we choose such a contour of integration which encloses both 
the nodes of the wave function G(r) and the boundary point, r = 0. Then the quantity N is 
depends on both the radial quantum number, n r , and the behaviour of the wave function near 
the origin, and is given by 



N = n r + yx 2 + W$ -V 2 , n r = n + {s + l)/2 , n = 0, 1, 2, ... (12) 

A further application of the theorem of residues to the explicit form of functions Ri(r) easily 
solves the problem of taking into account nodes of the wave functions for exited states. 

(2) Recursion formulae. We begin with investigation of behaviour of the functions Ri(r). 
From the system (6) we have 

i?o(r) = -^/m2-£ 2 , (13) 

where the minus sign is chosen from boundary condition. Then the function Ri(r) has a simple 
pole at the origin, owing to the Coulombic behaviour of the potentials at this point, while the 
function Rk(r) has a pole of the order k. Hence Rk(r) can be represented by the Laurent series 



oo 



R k (r)=r- k Y,Ry , (14) 

which makes it possible to write the quantization conditions (11) as 

R k k +1 = N5 k , . (15) 

Now, by analogy with Rk(r), let us also represent the functions Qk(r), involved in the 
expansion (5), as a power series in r: 

oo 

Qk{r)=r- 2 ~ k Y J Q> t ■ (16) 
From the equation Q(r)[E + m(r) — V(r)] = m'{r) — V'(r) we then obtain 



k-l 



Qi 



E + m 



3=0 3=0 



(17) 



Finally, substituting expansions (14) and (16) into the last equation of the system (6) and 
collecting coefficients of the like powers of r leads to the recursion relation in terms of the 
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Laurent coefficients, R^\ 



R! = 



2R» 



k-l 



k-2 i 

k-j \ ~> \ ~> r>3 ok-2-j 



k + mt 1 + E E W-7 -EE W-l 

j=l p=0 j=0 p=0 



~xQi 3 + $k,i E EjE k _j - 2E k ^ 1 V i - k+1 

i 

+S k ,2j2( V p V i-p ~ W P W i- P ) ~ S k,i 2mW t - 5 ifl 5 K2 X (x + 1) 

p=0 



(18) 



where for universality of designations we put = Rn, i?° = 0, i > 0. 

In the case i ^ k, this formula is intended for obtaining coefficients When % — k, by 
equating the explicit expression for R k +1 to the quantization condition (15) we get the recursion 
relation for the energy eigenvalues 



En 



m 



N 2 + Vr 



^ (NyjN* + V 2 - Wft 2 - Wo) ,k = 0, 



Eh — 



R l R l 



/k-l k 



i(EE R P R k k + -i 3 + 2e(A; -2)E K R U 

L^O \j=2p=0 



k 

E 

P =i 



k-2 
k 



2{E n R\ + V R° ) 

k—1 k 

-EE^G" J ! + s k ,i(v v 1 -w w 1 )-xQ K , 

j=o p =o ) (19) 

-4 (V 1 + E E K R U - E E Q 3 pRl-T j 

fc— 1 2 

+ £ - 2£ fc _ 1 y 1 + 4,2 E(^^- P - ^W 2 _ p ) 

j=l p=0 

-5 kA 2mW 1 - xQt 3 ^ ,k>0. 

Here E does be the exact solution to the Dirac-Coulomb equation [38,39] and we use the step 
function 

Q(k) =1, A;>0, 
= 0, k < . 

Thus, equations (18) and (19) determine the coefficients of perturbation expansions of en- 
ergy eigenvalues and eigenfunctions for screened Coulomb potentials in the same form both for 
ground and exited states. 



(3) Examples of application. As a check of the obtained formulae we calculate the energy 
eigenvalues for the pure-vector, screened Coulomb potential of general form. On applying the 
recursion relations (18) and (19), analytical expressions for the perturbation coefficients are 
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found to be equal to 



E 

E 2 
E 3 
E 4 



mN 

VN 2 + a 2 ' 
aVi , 

-^(3a 2 e- X ( X e + l)p 2 ) , 



r 



a\Ae 2 + 1) - a(2 X V + 3 X e + X ~ l)p 



(V 2 2 [a 4 e(5e 2 - 12) + a 2 e(Q X 2 p 2 - 5)p 2 



8p 6 
+X 
+V 4 

-3X(X 2 -I)(xe + 2) P 4 ]) , 



Ve(e 2 + 2) + X (4e 2 + 2) + 3e]p 4 
-5a 4 e(4e 2 + 3) + a 2 [Q X 2 e(2e 2 + 3) + 6 X (4e 2 + 1) - 25e]p 2 



8p 



■(v 2 V 3 [3a 4 {8e 4 - 20e 2 - 3) - a 2 [x 2 (32e 4 - 36e 2 - 10) 

+ X e(10e 2 - 24) + 9(6e 2 + l)]p 2 + X [30 X 2 e 3 
+24 X e 2 + lOe + X + X 3 (8e 4 + 8e 2 - l)]p 4 
+V 5 [-3a 4 (8e 4 + 12e 2 + 1) + a 2 [x 2 (16e 4 + 48e 2 + 6) 
+ 10xe(4e 2 + 3) - 15(6e 2 + l)]p 2 

+ [5 X 2 (4e 2 + 3) - 3 X 4 (4e 2 + 1) + 50 X e - 30 X 3 e - 12]p 4 ]) , 



(20) 



where a = Vq, e = E and p is \J\ — e 2 . 

One can verifies that the first three corrections coincide with those derived by McEnnan et 
al [24] with standard technique. 

The next example will be the attractive Yukawa potential, often utilized in relativistic cal- 
culations, which has not only the Lorentz-vector component, V(r) = — (a/r)e~ Xr , but the 
Lorentz-scalar term, W(r) = — {b/r)e~^ r , as well. Now the analytic expressions for perturba- 
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tion corrections to the bound state energy take the form 



Eq 

E 2 



m 



N 2 + a< 



(NVN 2 + a 2 -b 2 - ab) , 



a\ + bpe , 
1 



3a 3 e - a X (xe + l)p 2 + 2a 2 6(2e 2 + 1) 



E% — 



(21) 



4(6e + a)p 2 
+ab 2 e(e 2 + 2) 

+/i 2 [36 3 e 2 - b X ep 2 - b X 2 p 2 + 2b 2 ae(e 2 + 2) + ba 2 (2e 2 + 1) 

+3a 3 6e(2e 2 + 3) + a 2 6 2 (2e 4 + lie 2 + 2) - ab(3 X + (3% 2 - l)e)p 2 
+a6 3 e(3e 2 + 2)" 

+/i 3 [6 4 (-8e 4 + 13e 2 ) + 6 2 (3xe 3 + (2 X 2 + l)e 2 - 6*e - 5*V 
-36 3 ae(2e 4 - e 2 - 6) - 6 2 a 2 (4e 4 - 14e 2 - 5) 
-bae(3xe + 3x 2 - l)p 2 + &a 3 e(2e 2 + 3)" 
+X 2 p[9a%ep 2 - Qa 2 b 2 (2e 4 - e 2 - 1) 

+3a6 3 ep 2 (e 2 + 2) - 3ab X (x^ + l)p 4 



In order to assess the speed and accuracy of the perturbation technique for the Yukawa 
potential with various ratios of its component we consider energy eigenvalues for the pure 
vector case, V(r) = — (a/r)e~ Xr ; the pure scalar potential, W(r) = — (a/r)e -Ar ; and the 
equally mixed interaction, V(r) + W(r) = —(I + 7 )(a/2r) e~ Xr . Typical results of calculation 
are represented in Table 1 where the sequence of the sums of first terms from our expansion 
for relativistic binding energies is compared with the results, E num (in KeV), obtained by 
numerical integration. The calculation has been performed for s — 1, n r — 1, I — 1 and 
s = —1, n r — 1, I — states with parameters a = az, A = 1, 13az 1 ^, z = 74 ( a is the 
fine-structure constant and z the nuclear charge). 

As it can be seen from Table 1, in all cases we have two subsequences bounded below and 
above the energy eigenvalues. The average of these subsequences at the point of their maximal 
drawing together is proved to result in a quite good approximation to the exact value. 



To summarize, we have developed a semiclassical treatment of logarithmic perturbation the- 
ory for a Dirac particle in a central field. Based upon the ^-expansions and suitable quantization 
conditions, new handy recursion relations for solving the bound-state problem for the Dirac 
equation with the screened Coulomb potential having both vector and scalar component have 
been derived. Avoiding the disadvantages of the standard approach these formulae have the 
same simple form both for ground and exited states and provide, in principle, the calculation 
of the perturbation corrections up to an arbitrary order in the analytic or numerical form. And 
at last, this approach does not imply knowledge of the exact solution for zero approximation, 
which is obtained automatically. 
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Table 1 



k 




s — 1, n r 


= 1,1 = 1 






s 


= —l,n r 


= 1,1 = 






Ey 


Ew 


Ey+w 


Ey 


Ew 


Ey+w 


o 


20 


644616 


16.59173 


18 


292804 


20 


644616 


16.59173 


18 


292804 


1 


11 


091653 


7 348948 


10 


846326 


11 


091653 


7 348948 

1 • - 7 I\ 7 ' 7 -L V 7 


10 


846326 


2 


12 


415123 
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11 


810855 


12 


721342 


9.372375 


12 


003761 


3 


12 


264120 


8.784152 


11 


703753 


12 


500951 


8.995036 


11 


855249 


4 


12 


308677 


8 880691 


11 


730813 


12 


558797 


9 118034 


11 


890052 


5 


12 


292914 


8 838096 

v - • V 7 -.7 v7 V 7 '. 7 V 7 


11 


722213 


12 


537837 


9 063240 

'7 • V 7 V7 '.J ^ JL V / 


11 


878837 


6 


12 


299805 


8 860978 

v 7 • v7 V7 V 7 '.7 1 v7 


11 


725555 


12 


546910 


9 092488 


11 


883163 


7 


12 


296473 


8.847251 


11 


724110 


12 


542466 


9.074820 


11 


881275 


8 


12 


298236 


8.856234 


1 1 

1 1 


79/L7Q9 






9.086483 


11 


882175 


9 


12 


297242 


8.849965 


11 


724449 


12 


543482 


9.078246 


11 


881715 


10 


12 


297834 


8.854576 


11 


724631 


12 


544306 


9.084385 


11 


881963 


11 


12 


297465 


8.851033 


11 


724530 


12 


543784 


9.079601 


11 


881823 


12 


12 


297704 


8.853860 


11 


724588 


12 


544128 


9.083474 


11 


881905 


13 


12 


297544 


8.851529 


11 


724553 


12 


543893 


9.080233 


11 


881856 


14 


12 


297654 


8.853509 


11 


724575 


12 


544058 


9.083025 


11 


881886 


15 


12 


297576 


8.851782 


11 


724561 


12 


543939 


9.080555 


11 


881867 


E 

- L 'num 


12.297609 


8.852592 


11 


724567 


12 


543990 


9.081723 


11.881875 
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